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Abstract — In the present paper, a class of partial differential equations related 
to various plate and rod problems is studied by Lie transformation group meth- 
ods. A system of equations determining the generators of the admitted point Lie 
groups (symmetries) is derived. A general statement of the associated group- 
classification problem is given. A simple intrinsic relation is deduced allowing to 
recognize easily the variational symmetries among the "ordinary" symmetries of 
a self-adjoint equation of the class examined. Explicit formulae for the conserved 
currents of the corresponding (via Noether's theorem) conservation laws are sug- 
gested. Solutions of group-classification problems are presented for subclasses of 
equations of the foregoing type governing stability and vibration of plates, rods 
and fluid conveying pipes resting on variable elastic foundations and compressed 
by axial forces. The obtained group-classification results are used to derive con- 
servation laws and group-invariant solutions readily applicable in plate statics or 
rod dynamics. 



1. INTRODUCTION 

A wide variety of classical theories of plates and rods0 rest on linear fourth-order partial 
differential equations in one dependent and two independent variables. Some of them, such 
as the Poisson-Kirchhoff type theories for small bending of plates, are developed within the 
framework of the linear elastostatics to determine the state of equilibrium of thin elastic 
plates in terms of the transversal displacement of the plate middle plane, the derived gov- 
erning equations providing the background for solving problems concerning stability and vi- 
bration of such structural elements. Others (among them - the dynamic theory of Bernoulli- 
Euler beams, for instance) are deduced on the ground of the linear elastodynamics to describe 
the dynamic behaviour of rods in terms of the transversal displacement of the rod axis. 

The aim of the present paper is to study the invariance properties (symmetries) of the 
equations of the foregoing type with respect to local Lie groups of point transformations 
of the involved independent and dependent variables. The work is motivated both by the 
aforesaid wide applicability of the equations in question in structural mechanics, and by the 
remarkable efficiency demonstrated by the symmetry methods, especially when applied to 
differential equations arising in physics and engineering. 

Actually, once the invariance properties of a given differential equation are established, 
several important applications of its symmetries arise. First, it is possible to distinguish 
classes of solutions to this equation invariant under the transformations of symmetry groups 
admitted. The determination of such a group-invariant solution assumes solving a reduced 
equation involving less independent variables than the original one. Typical examples of 
group-invariant solutions are axisymmetric solutions, self-similar solutions, travelling waves, 
etc., which have proved to be quite useful in many branches of physics and engineering. 
For a self-adjoint differential equation another substantial application of its symmetries is 
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available. As is well known, the self-adjoint equations are the Euler-Lagrange equations of 
a certain action functional. If a one-parameter symmetry group of such an equation turned 
out to be its variational symmetry as well, that is a symmetry of the associated action 
functional, then Noether's theorem guarantees the existence of a conservation law for the 
smooth solutions of this equation. Needless to recall or discuss here the fundamental role 
that the conserved quantities and conservation laws (or the corresponding integral relations, 
i.e. the balance laws) have played in natural sciences, but it is worthy to point out that 
the available conservation laws (balance laws) should not be overlooked (as it is often done) 
in the examination of discontinuous solutions (acceleration waves, shock waves, etc.) or in 
the numerical analysis (when constructing finite difference schemes or verifying numerical 
results, for instance) of any system of differential equations of physical interest. It should 
be remarked also that the path- independent integrals (such as the well known J-, L- and 
M-integrals) related to the conservation laws are basic tools in fracture analysis of solids 
and structures. 

The aforementioned and many other applications of the symmetries of differential equa- 
tions and variational problems as well as the foundations of the Lie transformation group 
methods, including the basic notions, statements and techniques, can be found in the books 
by Ovsiannikov (1982), Ibragimov (1985), Bluman and Kumei (1989) and Olver (1993) (see 
also the references given in these books). In the present paper, however, as fare as the 
application of the symmetry groups of the equations studied is concerned, our attention is 
restricted to the constructing of group-invariant solutions and conservation laws. Of course, 
the first task is to find these symmetry groups, and as here we do not deal with a sin- 
gle differential equation but with a class of differential equations, this means to solve a 
group-classification problem. 

The layout of the paper is as follows. A detailed description of the differential equations 
to be studied as well as the variational statement for the self-adjoint equations among them 
are given in Section 2. Several examples of mechanical systems governed by such equations 
complete this Section. In Section 3, a system of equations determining the generators of the 
symmetry groups admitted by the equations of the class considered is derived, the invariancc 
properties inherent to the whole class because of its linearity and homogeneity being taken 
into account, and then the general statement of the associated group-classification problem 
is given. After that, the variational symmetries of the self-adjoint equations of the examined 
class are investigated. A simple intrinsic relation allowing to recognize easily the variational 
symmetries among the " ordinary" point Lie symmetries of such an equation is deduced, and 
explicit formulae for the conserved currents of the conservation laws corresponding to the 
variational symmetries via Noether's theorem are suggested. Group-classification results, 
conservation laws and group-invariant solutions arc presented in Section 4 for differential 
equations governing stability and vibration of plates of Poisson-Kirchhoff type. Similar 
results are displayed in Section 5 for equations governing vibration of rods on a variable 
elastic foundation and dynamic stability of fluid conveying pipes or rods compressed by 
axial forces. In the reminder of this Section, conservation laws for rods derived in the 
present contribution are compared to other ones reported in the literature. 



2. BASIC EQUATIONS 
Consider a fourth-order homogeneous linear partial differential equation 

A a ^ 5 {x)w aPl s + A a ^{x)w aPl + A afi {x)w aP + A a {x)w a + A(x)w = 0, (1) 

in two independent variables x = (x 1 ,^; 2 ) and one dependent variable w(x). Here and 
throughout: Greek indices have the range 1, 2, and the usual summation convention over 
a repeated index (one subscript and one superscript) is employed; w aia2 ... ak (k = 1,2,...) 
denote the fc-th order partial derivatives of the dependent variable, i.e. 

d k w 

w aia2 ... ak = dxaidxa2 _ dxa „ (* = 1,2,...)- 



Further, a similar notation will be used for the partial derivatives of any other function of the 
variables x l ,x 2 but, in this case, the indices indicating the differentiation will be preceded 
by a coma, e.g. 

a k Aaf3~/8 

AaPlS - ° (h. _ i o ^ 

The coefficients of equation ([!]) are supposed to be smooth functions possessing as many 
derivatives as may be required on a certain domain of interest, and to be symmetric under 
any permutation of their indices, i.e. 

Using the total derivative operators 

d d d d d 

D n = + Wc - h w n 



OX a OW OW^ OW^v OW^ua 

the equation (|l]) may be written in the form 

V[w] = 0, (2) 

where D is the linear differential operator given by the expression 

V = A a P jS D a DpD 7 Ds + A a(il D a D f) D 1 + A af3 D a D + A a D a + A. (3) 

An equation of form (|^) is the Euler-Lagrange equation associated with a certain varia- 
tional problem involving only one dependent variable if and only if the differential operator 
T> is self- adjoint, that is 

£> = £>*, (4) 
where T>* is the (formal) adjoint operator of T> (cf. Olver, 1993). The explicit form of T>* is 

V* = D a DpD 7 D s A al3 ' lS - D a Di 3 D 1 A af31 + D a D Si A afi - D a A a + A. (5) 

In such a case, W) can be associated with the variational problem for the functional 



A[w] = J \jV\w\dx 1 dx 2 



since the application of the Euler operator 



d d d d d 
3 D„- + D„D U - D^DvDv- + D^D l/ D (T D T - 
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on the Lagrangian density 



yields 



L=\wV[w], (6) 
V[w] = E(L), (7) 

due to (@) and (|). 

Let us give several examples of plate and rod structures whose governing equations are 
self-adjoint and belong to the class specified above. Henceforward, when regarding plates 
the variables x 1 ,x 2 will represent the coordinates of the plate middle plane. As for the rod 
problems, x 1 will be associated with the spatial variable along the rod axis, and x 2 - with 
the time. In both cases, the dependent variable w will represent the transversal displacement 
field. 

Example 1. Small bending of plates resting on elastic foundations. Consider a thin 
elastic plate of variable bending rigidity D(x) resting on an elastic foundation of Winkler 
type with variable modulus k{x) and subjected to an edge loading leading to the appearance 



of nonuniform membrane stresses N a @ (x) . In this physical situation, the equation governing 
the small bending of the plate assumes the following form 

A[DAw] - [(1 - v)e a ^e 0v D^ - N^]w^ + kw = , (8) 

the membrane stress tensor N a/3 being symmetric, N a @ = N@ a , and divergence free, i.e. 
N°!f — 0. Here: v is Poisson's ratio; A is the Laplace operator, that is A = 5 a ^d 2 /dx a dx^, 
where 5 a0 is the Kronecker delta symbol (<5 n = S 22 = 1, 5 12 = S 21 = 0) and e aP is the 
alternating symbol (e 11 = e 22 = 0, e 12 = —e 21 = 1). 

Example 2. Elastodynamics of Bernoulli- Euler beams. Consider a nonhomogencous 
Bernoulli-Euler beam with bending rigidity B(x 1 ) and inertia term H(x 1 ). The differential 
equation governing the small vibration of such a beam is (Chien at al., 1993): 

Bwim + 2i? iwm + B^iwn + Hw22 = 0. 

Example 3. Elastic beams resting on elastic foundations. Consider an elastic beam 
of constant bending rigidity K and constant mass density m (mass of the beam per unit 
length), resting on an elastic foundation with variable modulus k(x). Suppose it is subjected 
to a constant follower force p. Then, according to the Bernoulli-Euler theory, the differential 
equation for small transverse vibration of the beam is (see Smith and Herrmann, 1972): 

Kwim + pwn + k(x)w + mw22 = . (9) 

Example 4. Pipes conveying fluid. Consider an elastic circular-cylindrical pipe of 
uniform outer radius of the pipe cross section, which is supposed to be small in comparison 
with certain characteristic pipe length (for a simply supported pipe say the length of the 
span). Let the pipe conveys inviscid incompressible fluid with a flow velocity U = const. 
Then, the equation of motion is (see Gregory and Paidoussis, 1966): 

EJw nn + MU 2 w n + 2MUw 12 + (m + M)w 2 2 = , (10) 

where E is Young's modulus of the pipe material, J is the (axial) moment of inertia of the 
pipe cross section, m is the mass (constant) of the pipe per unit length, M is the mass (also 
constant) of the fluid per unit length. 

Combining and generalizing equations @ and ( |Io[ ) presented in Examples 3 and 4, in 
Section 5 we will pay particular attention to the differential equations of the form 

iwnxi + X a0 w aj 3 + k(x)w = , (11) 

where 7^0 and x°*^ are rea l constants, while k(x) is a smooth function. 



3. SYMMETRIES AND CONSERVATION LAWS 

Consider a local one-parameter Lie group of point transformations acting on some open 
subset 51 of the space R 3 representing the independent and dependent variables £ 1 ,x 2 ,u> 
involved in our basic equation (J2j) . The infinitesimal generator of such a group is a vector 
field X on R 3 , 

X = e{x,w)-^+i 1 {x,w)-^-, (12) 

whose components ^(x,w) and n(x,w) are supposed to be functions of class C°°on fL By 
virtue of Theorem 2.31 (Olver, 1993), a vector field X of form ([l2]) generates a point Lie 
symmetry group of equation (Q) if and only if there exists a function A depending on x, 
w and derivatives of w (that is a differential function) such that the following infinitesimal 
criterion of invariance, 

X (V[w]) - XV[w] = 0, (13) 

4 

holds; here X denote the fc-th prolongation of X (Ovsiannikov, 1982). 



The invariance criterion (|13J) leads, through the standard computational procedure (see, 
e.g. Ovsiannikov, 1982 or Olver, 1993), to the following results: 

(i) each equation of form (Q) being linear and homogeneous is invariant under the point 
Lie groups generated by the vector fields 

d d 
*o = w-, X u = u( X )-, (14) 

where u (x) is an arbitrary solution of the equation considered, the invariance criterion ( |l3|) 
being fulfilled with A = 1 for Xq, and A = for the generators X u ; 

(ii) an equation of form @ admits other vector fields (|l2|), in addition to the aforemen- 
tioned ( |l4| ) , if and only if they have the special form 

X = (15) 

the functions £ M (x) and a (x) being nontrivial solutions of the following system of determin- 
ing equations (called further the DE system for easy reference): 

pj^fhi + ( CT _ \)A afJ ~< s - A a ^^ - A af3 » s ^ - A a ^ 5 ^ - A^ 5 ^ = 0, (16) 

(17) 



4A a ^a tfl - 2A a(i ^ v C.^ - 2A a ^ v ^ llu - 2A^ V Q V 

+£»A°<Pi + {a- \)A af} ~< - A"! 3 ^^ - A a ^^ - A^Q = 0, 

6A^a^ v - 2A<*^°^ va - 2A^"°^ va 

+3A^«7, P - (3/2)A^^ u - (3/2)A^X»» ( 18 ) 
+ ^Af + (a - X)A a/3 - A a »f - A^£_ a = 0, 



a Aafj.ua _ AfivaTta 



-3A a ^ V 0- au - A^ va C 



2A<*»v tl > - A^Q V 



(19) 



fJLV 

+(*A% + \)A° - = 0, 



A^a^s + A a ^u, a01 + A af} a^ & + A a a, a + + (a - X)A = 0, (20) 

for a certain function A depending on x 1 and x 2 only. (Here, by a trivial solution we mean 
not only £ M — 0, a — 0, but also £ p = 0, a = c = const ^ 0, since the latter leads to the 
vector field cXq generating the same group as Xq which is already identified to be admitted 
by each equation of the type considered.) 

Thus, given an equation of form (0), the question is whether there exist vector fields 
X 7^ cXq of form ( jllf ) which leave it invariant, and the answer depends on whether the 
respective DE system has at least one nontrivial solution. In this context the coefficients 



of the equation are supposed to be known functions, and thereby (16) - (J20|) constitute an 
over-determined system of linear homogeneous partial differential equations with respect to 
the unknowns £^ and a. Therefore, as a rule, it turns out possible to find in an explicit form 
some (or even all) nontrivial solutions of the DE system, and thus to determine several (all) 
additional point Lie symmetry groups inherent to the equation in question. 

It should be remarked that various equations of form ^) admit only the point Lie groups 
generated by the vector fields ([l4]) with u(x) being any solution of the respective equation. 
For instance, it is easy to check that all equations of the form (O) such that x Q/3 = S al3 
and k(x) = p(x), where p{x) is an arbitrary polynomial of x 1 and ar, belong to this variety. 
Without too much difficulties one can ascertain that the same holds true for the equations 
of the form (|) with D = const, N af3 = 5 a0 and k(x) = p{x). 

On the other hand, there are equations of the foregoing type which are invariant under 
a larger group; an immediate example is the biharmonic equation, A 2 w — 0, which admits 



the seven-parameter group generated by the linear combinations of Xq and the following six 
additional basic vector fields (cf. Ovsiannikov, 1972): 



d v d , d v _ n i 2 a 



dx 1 ' ' 9a; 1 9a; 2 9a; 

9 „ n 9 n 9 



*i = £t, Xa^x^-a; 1 ^, X 5 = 2zV-^ - (x 1 ) 2 - (* 2 ) 2 



9 _ , 9 



A 2 = -^, X 4 = a; 1 -^ T +a; 2 -^2, A 6 = (a; 1 ) 2 - (a; 2 ) 2 + 2a; 1 a; 2 ^2 + 2a; 1 w^. 
9ar ox L 9ar L v J da; 1 9ar Ow 

An important problem naturally arises in the light of the above note. It may be placed 
in the category of the so-called group-classification problems (see Ovsiannikov, 1982) and 
consist in determination of all those equations of the type considered that admit a larger 
group together with this group itself. Its most general statement assumes all functions 
A a Pv s (x), A a ^{x),A al3 (x), A a (x), A(x),£ a (x) and a(x) involved in the determining equa- 
tions (|l^) - (|2^) to be regarded as unknown variables and to find all solutions of this system. 
Here we are not going to study this rather complicated nonlinear problem in general. How- 
ever, in Sections 4 and 5, restricting our attention to the equations of form (||), D = const, 
and (|ll|), respectively, we will examine the corresponding group-classification problems. 
Let us now specialize to the case of self-adjoint equations of form (0). Suppose that 

V[w]=0, T> = V*, (21) 

is such an equation. Then, of particular interest are its variational symmetries - the Lie 
groups generated by the so-called infinitesimal divergence symmetries (see Definition 4.33 
in Olver, 1993) of any variational functional with ( |2~l| ) as the associated Euler-Lagrange 
equation. (Note that if two functionals lead to the same Euler-Lagrange equation, then they 
have the same collection of infinitesimal divergence symmetries.) This interest is motivated 
by the fact that, in virtue of Noether's theorem, each variational symmetry of a given self- 
adjoint equation corresponds to a conservation law admitted by the smooth solutions of the 
equation. Thus, if a vector field X of form ( |l2| ) is found to generate a variational symmetry 
of equation (pl|), then Noether's theorem implies the existence of a conserved current, which, 
in the present case, is a couple of differential functions P a such that 

D a P a = QV[w], (22) 

where Q is the characteristic of X; by definition 

Q = tj (23) 

The total divergence of the conserved current P a vanishes on the smooth solutions of ( pl| ) 
and so we have the conservation law 

D a P a = 0, (24) 

22|) being its expression in characteristic form, and Q - its characteristics. Therefore, to 
derive the conservation laws of the foregoing type, one can proceed by first determining the 



i ■ 1 2x W 7T> 
Ox Ow 

9 „ i 9 9 „ i 9 



variational symmetries of equation (21), and than using their characteristics (23) to find, 
from (p2]), explicit expressions for the corresponding conserved currents. 

Having analyzed earlier the invariance properties of the whole class of equations (EH , it 
is convenient to base the determination of the variational symmetries of equation (|2lf ) on 
the following observation. A vector field X of form (|l2| ) generates a variational symmetry 
of equation ( pl| ) if and only if X is an infinitesimal symmetry of this equation, that is ( |l3[) 
holds, and 

dn 
dw 



x (PM) +\^r.+ D ^ ) p N = o. (25) 



This is a consequence of Lemma 4.34 and Proposition 5.55 (Olver, 1993), see also Lemma 
7.46 (Olver, 1995). Subtracting (|l3|) from p3) we can replace the latter with 



^L+D^e + X) V[<r] =0. 



and as D[w] is not supposed to vanish identically we arrive at the conclusion that 



|U^ + A = 0, (26) 

is a necessary and sufficient condition for an infinitesimal symmetry admitted by a self- 
adjoint equation of form (^) to be its infinitesimal variational symmetry as well. It should 
be remarked that the same holds true for any self-adjoint partial differential equation in 
one dependent variable w and n independent variables x — (x 1 , . . . ,x n ); of course, in the 
general case the summation index fi will take the values 1, . . . , n in both formulae (|12| ) and 
(p6|). For a vector field of form ( jilt ) the relation ( p6|) simplifies, and reads 

<r + f£ + A = 0. (27) 

Thus to find the variational symmetries of an equation of form (^l|) , it suffices to check 
which of its "ordinary" symmetries satisfy the additional requirement (p6|). For instance, 
the result (i) implies that Xq = wd/dw does not generate a variational symmetry of any 
equation of form (pi]), while a vector field X u = u(x)d/dw generates a variational symmetry 
of an equation of form ( ]2l| ) whenever u(x) is its solution (this is a common property of all 
systems of linear homogeneous partial differential equations, see Section 5.3 in Olver, 1993). 

Suppose one has established that a vector field X with characteristic Q generates a vari- 
ational symmetry of a given equation of form and now wishes to find the conserved 
current P a of the corresponding conservation law (p4|). For this purpose one can use for- 
mulae (5.150) and (5.151) given by Olver (1993) which express (in an explicit form) a null 
Lagrangian as a divergence. Indeed, in this case the right hand side of ( p2| ) is a total diver- 
gence or, in other words, a null Lagrangian. However, bearing in mind the recommendation 
of Olver (1993) to use these formulae only as a last resort since "the homotopy formula 
(5.151) can rapidly become unmanageable", in the present paper we suggest another way 
for determination of the sought conserved currents. 

Our starting point is the so-called Noether identity (cf. Ibragimov, 1985): 

X OC) + (D a C)C = QE(C) + D a N a (£), (28) 

oo 

which holds for any differential function C and vector field X of the types considered here. 
In (|28|), N a are the differential operators given by the expressions 



r>1 [ UUJ an 1 ...n r s>1 uuJotn 1 ...n r v l ...v tl 

Q = r) — (, a w a being the characteristic of the vector field X. Setting C = L in (|28|), and 
taking into account (|^) and (R), after a little manipulation we obtain the identity 

D^N^{-wV[w]) = -wX {V[w}) -{f]+ (D^)w - 2Q}V[w], (30) 

4 

valid for any self-adjoint differential operator V of form @ and vector field of form (|l2|). 
In particular, for X v = v(x)d/dw, where v{x) is an arbitrary smooth function, we have 

^ = 0, Q = V = v, (31) 

and hence 

X v [V[w]) = V[v], (32) 

4 



since T> is a linear differential operator. Substituting (|3l]) and ( p2[ ) into ( |30| ) we obtain 



D^N^i-wVlw]) = vV[w] - wV[v], (33) 

which is nothing but the reciprocity relation associated with the equation D[w] = 0. Under 
the additional assumption v — u(x), where u(x) is an arbitrary smooth solution of the latter 
equation, the reciprocity relation d33j) becomes 



D^i-wVlw]) = uD[w\. (34) 

Taking into account (|34|) we can give now the following general formula for the conserved 
currents P a of the conservation laws with characteristics Q = u corresponding to the in- 
finitesimal variational symmetries X u = ud/dw of equation (fjl[): 

P a = P" u) + G a , 

where 

P ( «) = N a (-wD[w]), (35) 
and G a is any null divergence. Of course, 

D^Pfa = uD[w], 

and 

D»P( U) = 0, (36) 

on the smooth solutions of the equation (pl|). 

Next, let X be an infinitesimal variational symmetry of equation ( ^l|) with characteristic 
Q = wo — w^^ 1 . Then, on account of (^5|), ( ^o|) takes the form 

f-^wV[w]j = QV[w], 

and hence we can write down the following explicit formula for the conserved currents 
P a of the conservation laws with characteristics Q = wo — w^ corresponding to the 
aforementioned variational symmetries of equation (21), namely 

P a = B a + G a , 

B a = N a (-^wD[w}) + (wCA^DpD^sw - w^A^D^Dgw) , (37) 
where, as before, G a is any null divergence. Of course, 

= QV[w}, 

and on the smooth solutions of respective equation ( pl| ) we have 

= 0, 

Let us remark that the special null divergence 

^ (wCA^DpD^Dsw - weA^DpD^Dsw) , 

is used in the expression fl37|) for the conserved current B a to cut the fourth-order derivatives 
of the dependent variable w away since in practice we are usually interested in conserved 
currents which involve derivatives of order not higher than k — 1, where k is the order of the 
equation considered. Making use of (|9|) it is easy to check that the right-hand side of (37) 
incorporates derivatives of w of order less than fourth. In the subsequent Sections just (35) 
and ([37]) will be referred to as the expressions for the conserved currents of the conservation 
laws with characteristics Q = u and Q = wo — respectively, derived for equations of 

the form (12 



To summarize, given an equation of form (|2l]), the crucial point on the way of deriving 
conservation laws admitted by its smooth solutions is to find vector fields of form ( |l5| ) 
generating "ordinary" point Lie symmetries of the given equation. For that purpose, we 
should look for solutions of the respective DE system ( |l~6| ) - ( p0| ) . Once such vector fields 
are found, it is easy to check which of their linear combinations satisfy the requirement 
( p6[ ) and hence generate variational symmetries of the equation considered. Now, using the 
characteristics of these symmetries we first construct the operators N a from formulae ( p9[ ) 
and then calculate from (37) the conserved currents of the corresponding conservation laws. 



4. SYMMETRIES, CONSERVATION LAWS AND GROUP-INVARIANT SOLUTIONS 

OF PLATE EQUATIONS 

In Section 2 (Example 1), we have quoted the self-adjoint equation (||) describing the 
small bending of a plate resting on an elastic foundation. Many problems concerning stability 
and vibration of isotropic thin elastic plates are studied on the ground of this type of 
equations. Here, we analyze the invariance properties of a generic equation of this form 
under the assumption that the bending rigidity of the plates considered is uniform, that is 
D = const. In this case (g) may be written as follows 

A a ^ s w aPjS + A afi {x)w aP + A(x)w = 0, (38) 

with 

Aat 3-yS = l^gafipS + g^gfJS + gctifiPr^ A J = Q) (39) 

3 

assuming that 

A a P = l-N af3 , A=—k. 
D D 

In view of the general results of Section 3, it is clear that X u — u(x)d/dw generates 
a variational symmetry of any equation of form ( |38|) whenever u{x) is its solution, while 
Xq = wd/dw alone could never generate a variational symmetry of an equation of form 
(j38[), though it always is its infinitesimal point Lie symmetry. Substituting ([59]) into the 
determining equations (|l6|) - (|20| ) and taking into account that A a ^ = 0, A a — 0, we obtain, 
after a straightforward computation, that an equation of form ( |38|) is invariant under a point 
Lie group generated by a vector field X of form (|l5|), X ^ cXq, if and only if 



1 



2 C> (4°) 

<^ef M + ^C-^C = °. (4i) 

eAf - A a ^ - A^Q + 2^A a P + 2S aT S^^ TV = 0, (42) 

A av ^ v - A^i% = 0, (43) 

2eA IM +A^A + A^Xr„u = 0- (44) 

At that, 



A = -^C- (45) 



3 

Substituting expression ( |40| ) into (|l^) , and expressions @) and (pt5| ) into condition (^) , we 
immediately arrive at the conclusion that the generator of such a group is a vector field of 
form 

each such symmetry of ( |38| ) being variational symmetry of the latter equation as well, and 
hence there exist a conservation law with characteristic Q — {1/2)^ w — and conserved 
current B a given by (p^) admitted by the smooth solutions of the equation considered. 



Thus to derive the conservation laws, which correspond to the variational symmetries of an 
equation of form (|3|) it suffices to know the results of the group classification of the class 
of equations in question; of course, the same holds true for the derivation of group-invariant 



solutions to (38). This group-classification problem is studied in Vassilev (1988), (1991) and 
(1997). The classification results presented below are obtained in these works. 



It is shown that the scalar fields 



A" v 6„ 



1/2 



'(2) 



s (3) 



(1),M S (1) 



1/3 



(47) 



are of key importance for the group classification of the considered class of equations. These 
scalar fields are called the invariants of equation ([38j) since here they play a role similar to 
the role that Laplace's and Cotton's invariants play in the group classification of the second- 
order linear partial differential equations (see Ovsiannikov, 1982 and Ibragimov, 1985). The 
following two properties of the scalar fields (E^) give us both an additional reason to call 
them invariants of ( |38| ) and explicit expressions for the invariants of groups admitted by 
First, if an equation of form (|38|) admits a vector field of form (E6I), then 



= 1,2,3), 



and hence Uu\ — 



w \J s {j) are invariants of the corresponding Lie group whenever S(j} ^ 0. 
Second, if an equation of form (^8|) admits a vector field of form ( |46| ) and is such that at least 
two of its invariants (^7|), say and sm (k ^ I < 3) , are not identically equal to zero, then 
s (k)/ s (i) is an invariant of the corresponding symmetry group. Note, that the invariants s/^n 
and sm of such an equation of form ( ^8|) provide two couples of functionally independent 



invariants, namely U^) 



and S(fe)/s(() as well as Un\ 



and 



(fc)/ 6 '(0> 



of 



the admitted symmetry group, both couples being readily applicable for constructing group- 
invariant solutions to the respective equation. However, if even one of the invariants ( [47| ) 
of an equation of form ( |38| ) is not identically equal to zero, then this equation admits at 
most a 3-parameter group with generators of form (Q) . On the other hand, if all invariants 
(E^) of an equation of form ( |38| ) are identically equal to zero, then this equation admits a 
6-parameter group with generators of form (|6j). Below, the latter case is set out in detail. 

Let w(z) ^ const be an analytic function of the complex variable z = x 1 + ix 2 , and let 
E u be the equation of the form (Ssh with coefficients 

A 11 = -A 22 = 4Re{0} , A 12 = A 21 = - 

where d> is the Schwarzian derivative of the function uj, 



4Im{0} , 
that is 



.4 



(48) 



(49) 



<f> is the complex conjugated of <f>, and the prime is used to denote differentiation with 
respect to the variable z. Substituting ( ^§1) into (|47| ) one can see that all invariants of E w 
are identically equal to zero. Then, taking into account the DE system (^l|)-(^4|), (|4S| ) and 
(^9|), one can verify by direct computing that E u admits the 6-parameter group generated 
by the vector fields 



d 



1 







the functions £^ 



^dx^ 2 S ^< dw 
being given by the expressions 



-,6), 



I 13 


= Re{wi} , 


P 2 


= Im{wi} , 




= Re {iwi} , 


P 2 

5(2) 


= Im{iwi} , 


£(3) 


= Re{w 2 } , 


£(3) 


= Im{w 2 } , 


C(4) 


= Re{iLU 2 } , 


P 2 
5(4) 


= Im {iuj 2 } , 


£(5) 


= Re{w 3 } , 


P 2 
S(5) 


= Im{w 3 } , 


C(6) 


= Re{iw 3 } , 


P 2 
4(6) 


= lm{iu 3 } . 



where 



wi = — , uj 2 = — W 3 = — . 

a/ u/ a; 



(50) 



It should be remarked that each equation of form ( |38| ) which admits a 6-parameter group 
with generators of form ([l6|) is of type E u , meaning that it can be generated in the above 
manner using a suitable analytic function w. The coefficients of each equation of this type 
are of the form A af} = 5 a ^ v ip iliv , A = (l/8)S ati 5^ip ta0 (p^ u , where ip is a harmonic 
function, that is 5 a ^(p Qj3 = 0, and vice versa. It is noteworthy that each equation with 
variable coefficients of type can be mapped to an equation with constant coefficients 
belonging to the same family. It is easy to verify by direct computing that the equation 
Eu corresponding to an analytic function u> whose Schwarzian derivative is not constant 
transforms to a constant coefficients one under the following change of the independent and 
dependent variables: 

y a =f a (x\x 2 ), W = wU(x 1 ,x 2 ), 



/- X £feL U( X \x 2 ) 



f 1 (x\x 2 ) =Re|y f^dzj, f 2 (x\x 2 ) =Im 
where / is any linear combination of the functions (|5Cj) such that / ^ 0, i.e. 

/ = kiLUi + k 2 LU 2 + k 3 uj 3 , 



ff 



(51) 

1/2 



(52) 



where k±, k 2 , and k 3 are complex constants such that k\ + k\ + k 2 ^ 0. 

Consider, as a simple example, the equation E u corresponding to u = z^^^ 2 , where x 
is a positive real constant. In this case (^) gives 4> — 4 (2 — x) z 2 and hence, according to 
i8h, the coefficients of En, read 



A 11 = -A 22 = (2 - x) 



(x^-(x 2 ) 2 



, A 12 = A 21 = (2 - x) - 



.21 



(xi) 2 + (x 2 ) 2 



2x 1 x 2 



(xi) 2 + (x 2 ) 2 



A = (2 - x) 



Or 1 ) 2 + (x 2 ) 2 



Using the function / = z obtained from J52] ) for cj = zV^/ 2 , fc 1 = fc3 = 0, k 2 = 1 
we introduce, according to (51), the new independent and dependent variables 



1 



In 



(a 1 ) 3 + (*»)* 



y = arctan 



W = w 



(^) 2 + (x 2 y 



-1/2 



(53) 
(54) 



Note that the inverse transformations are given by the expressions 



x 1 = e y cosy 2 , 



2 7/ • 2 

x = e" sm y , 



w = W 



1/2 



(55) 



Under the change of the variables according to (54), the considered equation E w transforms 
to the following one, 



d 4 W 



d 2 W d 2 W 
— X ~\- X 

dy a dy l3 dy fl dy l/ dy 1 dy 1 dy 2 dy 2 



^>c 2 W = 0, 



(56) 



which belongs to the same class (since it corresponds to the analytic function lu = e z V^/ 2 ) 
but whose coefficients are constant. 



Equation (|5E|) admits the 6-parameter group of variational symmetries generated by the 
basic vector fields 



Va 

v 3 

V 5 
V 6 



d 

dy a ' 

«» (Oy 2 ) A + sin (Oy 2 ) A + e ^cos (% 2 ) A, 



-e e ^ cos (% 2 ) A + e 9 " 1 cos (% 2 ) A _ Qe^wtm (Oy 2 ) A, 
- 9yl «»B (% 2 ) A _ e -tf sin JL _ fe-^^oo. (^) |_, 



^ cob (% 2 ) A + e -V sin JL _ ^ ttfl01 (£tf») A, 



<9w ' 
_d 

dw ' 

1 



where = y/jt/2. These vector fields give rise to six linearly independent conservation laws 
for equation (p6| ). The characteristics of these conservation laws are 

= 5 w ^r y J W ~ W M ^ 0' = i, • • • . 6) . 

Here, Vj are regarded as operators acting on the functions £ : i? 2 — * i?. The corresponding 
conserved currents can be easily calculated from (37). 

Finally, let us remark that each one-parameter group generated by a linear combination 
of the basic vector fields Vj can be used for constructing group- invariant solutions of equation 
(p6[). Consider, for instance, the group H (V3 + V5) generated by the vector field V3 + V5. 
The functions s = sin (Oy 2 } / cosh (Oy 1 } and u — WJ cosh (By 1 ) constitute a complete set of 
invariants for this group and hence, following the well known algorithm (Ovsiannikov, 1982; 
Olver, 1993), we seek the H (V3 + V^-invariant solutions of equation (|5^) in the form 



W = u (s) cosh (Oy 1 ) , s 



cosh (Oy 1 ) 

Substituting ( p7| ) into (p>6|), we get the reduced equation 
The general solution to this ordinary differential equation is 



(57) 



u (s) = Ci + C 2 In 



5+1 
S- 1 



C3S -I- C4S In 



8 + 1 

s- 1 



where Ci, C2, C3 and C4 are real constants. Hence the if (V3 + V5)— invariant solutions of 
equation (pq) are given by the expression 



W (y\y 2 ) = Ci cosh (fly 1 ) + C 2 cosh (fly 1 ) In 



sin 



{Oy 2 



cosh I 



C3 sin 1 



■C A 



sin 1 



cosh (Oy 1 ) 



In 



sin (fly 2 ) — cosh (Oy 1 

sin (fly 2 ) + cosh (Oy 1 ) 
sin (fly 2 ) — cosh (Oy 1 ) 



Using the inverse transformations ( p5|) one can convert the above solutions of equation ( p6| 
into solutions of the equation E^, ui — zV^, with variable coefficients (|53|). 



4. SYMMETRIES, CONSERVATION LAWS AND GROUP-INVARIANT SOLUTIONS 

OF ROD EQUATIONS 



In Section 2, combining and generalizing Examples 3 and 4 we have introduced the class 
of self-adjoint partial differential equations 



with coefficients 



A 1111 w 1111 +A al3 w a!3 + Aw = 0, (58) 

A 1111 = 7, A^= x a \ A = k(x), (59) 

where 7 = const ^ 0, x Q/3 are arbitrary constants (but (x 12 ) 2 + (x 22 ) 2 7^ 0j otherwise ( ^8|) 
degenerates and becomes an ordinary differential equation), and k(x) is an arbitrary func- 
tion. Equations of this special type are used by many authors to study applied engineering 
problems concerning dynamics and stability of both elastic beams resting on elastic foun- 
dations (see e.g. Smith and Herrmann, 1972) and pipes conveying fluid (see e.g. Gregory 
and Paidoussis, 1966). In the present Section we first examine the point Lie symmetries of 
( p8[ ) and solve the corresponding group-classification problem. Then we derive conservation 
laws and group-invariant solutions of various rod equations of form (j58[). 

Consider the group-classification problem. In view of the results (i) and (ii) of Section 3, 
it is clear that each equation of form ( |58| ) is invariant under the point Lie groups generated 
by the vector fields Xq — wd/dw and X u = u (x) d/dw where u (x) is any smooth solution 
of the respective equation and the objective is to find those equations of the type considered 
which admit vector fields X of form ([l5]), X 7^ cAo, c = const ^= 0. 

Substituting @ into (|l6|) - (|2fj), and taking into account that A a , A a ^ and all A a ^ s 
except A 1111 are equal to zero, we obtain after a little manipulation the following system of 
determining equations for the functions £ M (x) and a (x) associated with the sought vector 
fields X of form @: 

C 2 ! = 0, 2a A - 3£n = 0, (60) 

5<m + 2x"^i - 2x 12 ^ 2 = 0, (61) 

X 22 (2^-e 2 2 ) = 0, X 12 (3e:!-e 2 2 )-X 22 e^0, (62) 

2x 12 ^,2 - X 22 ^ 22 = 0, 2 X 12 a,! + x 22 (2a, 2 - £ 2 22 ) = 0, (63) 

aim + X Q Va/3 + + = 0, (64) 

the auxiliary function A being expressed by 

A = a-4£\. (65) 

We look for the equations of form (^8|) whose coefficients 7, x"^ an d K (x) are such that 
system (Eol) - (|64|) possesses solutions different from the trivial one = 0, a = const =/= 0. 



Observing the system of determining equations (60) - ( p4j ) we see that for the purposes 
of the group-classification it is convenient to divide the equations of form ( p8| ) into three 
subclasses depending on whether the coefficients °f the given equation are such that: 
(A) x 22 ^ 0,det(x a/3 ) ^ 0; (B) x 22 ^ 0,det( X Q/3 ) = or (C) x 22 = 0,det(x Q/3 ) ^ 0. This 
covers all possibilities except for x 22 = 0,det(x Q/3 ) = when ( ^8|) becomes an ordinary 
differential equation that falls outside our interest in the present paper. 

For convenience we introduce the notation 



d v _ ( x , X 12 2 \ 9 2 d _ ( , x 11 2 \ 9 2 8 
ra" V x y ox L ox z \ x / 9x L ox z 



Let x 7^ 0. Then, the first equation (p2J) is equivalent to 



2^i " C 2 2 = 0. 



(66) 



Differentiating ( |66[ ) with respect to x 1 and taking into account the first equation (60) we 
obtain £ = 0. Hence, ( |6~L| ) and the second equation (62) reduce to 

X i %~X 1 % = 0, x 12 t\ - X 2 % = 0. (67) 

Consider the subclass (A): \ 22 ^ 0, det(x Q/3 ) ^ 0. Then, the first equation ( |60| ) together 
with ( |66| ) and ( |67| ) lead to ^\ = £\ = ^\ = £ 2 2 = 0, i.e., £ Q = c" = const. Consequently, 
the second equations of (J6TJ) and ( p3| ) imply a = c = const. At this, the first equation ( |63| ) 
is satisfied and (|64|) becomes 

c a n, a = 0. (68) 

All this means that when \ 22 an d det(x Q/3 ) ^ the DE system has only the trivial 
solution unless 

k(x) = /(/3V -(3 l x 2 ), (69) 

for a certain smooth function / ^ const and certain constants (3 a such that + 7^ 0. 
In this latter case, the DE system has the nontrivial solution 

e=p\ e = P 2 , v=o, (70) 

and so the differential equations of that kind admit additionally the one-parameter symmetry 
group associated with the vector field /3 Y\ + fi 2 Yi- In the case K — const, ( |68| ) is satisfied 
for any couple of constants c Q , and hence such equations admit the 2-parameter symmetry 
group with generators Yi and Y2. This completes the analysis of subclass (A). 

Consider now the subclass (B): \ 22 7^ 0,det(x"' 3 ) = 0. Differentiating the second equa- 
tion ( |67j ) successively with respect to x 1 and x 2 we obtain 

X X S,12 — u > X 5,12 X <i,22 — u ; 

which, on account of £ xl = 0, leads to £ 12 = £ 2 2 = 0- This result, together with (^) imply 
£ 2 22 = 0, and the nontrivial solution of (pOj) - (|63|) is now obvious: 



e=c l +c*(x l + ^x^j, e = c 2 + 2c 3 x 2 , a = 0, (71) 

with c 1 , c 2 and c 3 - arbitrary constants. Equation (|64|) reduces to 

^ + ^ = 0. (72) 

For an arbitrary k(x) it leads to = 0; it is easily verified that for k(x) = x 1 +x 2 x 2 , £ Q = 
is the only solution of ([72]). 

If k(x) is a function of form (|6TJ), then ( |72| ) implies ( |70| ) as a nontrivial solution to the 
system ( |60| ) - (|64|). Therefore, a generic equation of this kind admits only the one-parameter 
group generated by f3 1 Yi + f3 2 Y2, unless k(x) is of one of the following two special forms. 
The first one is 

k(x) — kq (f3 + a; 2 ) 2 , kq — const ^ 0, j3 = const, (73) 

when the nontrivial solution of (^Tj|) - (Q) is (|7l]) with c 2 = 2/3, c 3 = 1, c 1 - arbitrary, and 
hence the equations of subclass (B) with n(x) of form (|7^) admit the 2-parameter symmetry 
group generated by the vector fields Y\ and 2/3Y" 2 + Y 3 . The second special form of the 
function n(x) is 

/ x i 2 . \~ 4 
k(x) = ko ( /3 + 22 x2 ) ' K ° = cons t 7^ 0, /3 = const, (74) 

when the respective differential equations admit the 2-parameter group spanned over the 
vector fields f3Y x + Y 3 and (x 12 /x 22 ) Y t +Y 2 . 



Another extension of the symmetry group is possible if there exist two constants /3 1 and 
/3 2 , as well as a smooth function / ^ such that 

K{x) = (/3 2 + x 2 ) 2 f(y), y=(p 2 + x 2 y 1/2 (> + - ^x 2 ^j . (75) 

If k(x) is of form ( |75| ) , then the nontrivial solution of the determining equations ( |60| ) - ( |64|) 
is @ with c 1 = /3 + 2 (x 12 /x 22 ) /3 2 , c 2 = 2/3 2 , c 3 = 1 and the differential equations of that 
sort admit the one-parameter group, generated by 

/3 1 + 2 ^22 ^) Y l + 2 / ?2y 2 + *3 

only, except for the cases f(y) = hqV 4 ( k o — const), when k(x) takes the form (]74|), and 
f(y) — const ^ when k(x) becomes (J73| ) . 

The differential equations of form (pq) with x 22 ^ 0, det(x Q ' 3 ) = admit the 2-parameter 
group generated by Y% and Y% when k(x) = const ^ or the 3-parameter group with 
generators Yi, I2 and I3 when k(x) = 0. 

Finally, consider the subclass (C): x 22 = 0,det(x Q ' 3 ) ^ 0. Substituting x 22 = in the 
determining equations ([}0|) - (|63|), the latter simplify to 

e 2 i = 0, £ u = 0, x"6 - X 12 ^ = 0, 3^~e 2 2 -0, ^ = 0, ct i2 = 0, 
and their nontrivial solution is easily obtained: 

e=c 1 + c*L 1 +'^x 2 \ e = c 2 + 3c 3 x 2 , a = 0, (76) 

where c % are arbitrary constants (note that if x 22 = 0, then det(x"^) ^ assumes x 12 7^ 0). 
Equation (H) takes the form ( |72| ) and for an arbitrary k(x) it leads to £ Q = 0. If k(x) is 
of the form (|69|), such equations admit the one-parameter group generated by (3 l Yi + /3 2 Y2 
only, except for two special forms of k(x), namely 

k(x) = k ([3 + x 2 )~ 4/ \ (77) 

and 

k(x) = k ({3 + 2x l - . (78) 

where kq and f3 are constants. In the case ([7?]), the nontrivial solution of the determining 
equation (^) - ( p4] ) is (f76|) with c 2 = 3/3, c 3 = 1, c 1 - arbitrary, and the differential equation 
considered admits the 2-parameter group spanned over the vector fields Y\ and 3/312 + Ki- 
ln the case (p§) the group admitted is also a 2-parameter one, but generated by f3Y\ + I4 

and (x 11 /x 1 T y i+ 2y 2- 

Another extension of the symmetry group is possible if 

<x) = (/3 2 + a; 2 )" 473 f(y), y=(p 2 + x 2 y 1/3 (> + 2x' ^x 2 ^J , (79) 



where /3 Q are constants and / ^ is a smooth function. In this case, the nontrivial solution 
of © - (g) is © with c 1 = /3 1 + 3(x n /x 12 )/3 2 , c 2 = 6/3 2 , c 3 = 1 and we conclude that 
the equations of subclass (C) with k(x) of form (|79| ) admit only the one-parameter group 
generated by [/3 1 + 3(x u /v")j9 2 ]Yi + 6/3 2 F 2 + F 4 , except for /(y) = k 2T 4 (ko = const), 
when ( [79| ) coincides with J78| ) or /(?/) = const 7^ when k(x) has the form (|77|). 

Evidently, the differential equations of form ( p8| ) with x 22 = 0, det(x"' 3 ) 7^ admit: the 
2-parameter group with generators Y\ and Yi when k(x) — const ^ 0, and the 3-parameter 
group associated with Y\, Y 2 and Y4 - when k(x) = 0. 

The results of the above group-classification analysis are summarized in Table 1, where 
the equations invariant under larger groups are given through their coefficients together with 
the generators of the associated symmetry groups. 



Table 1. Equations of form (pa) invariant under larger symmetry groups. 



41 
ft 




npnprfltnrQ 


1 


k(x) = f(f3 2 x 1 -/3 1 x 2 ) 


/S 1 ^ + (3 2 Y 2 


2 


X 22 ± 0, det( X a ?) = 0, k(x) = (/? 2 + x 2 y 2 f(y), 
y = {/3 2 + x 2 y U2 [/3 1 + x 1 - (x 12 /x 22 ) x 2 } 


[f3 1 + 2(x 12 /x 22 )P 2 }Yi 
+2f3 2 Y 2 + Y 3 


3 


X 22 = 0, det(x^) ± 0, k(x) = (/3 2 + xY 4/3 /(»), 
y = (/3 2 + .x 2 )" 1/3 [/3 1 + 2.x 1 - ( X n /X 12 ) A 


IP 1 + 3(x n /x 12 )/3 2 ]n 
+6(3 2 Y 2 + 2F 4 


4 


X 22 + 0, det(x Q/3 ) = 0, k(x) = k (j9 + a; 2 )" 2 , 


Y u 2f3Y 2 + Y 3 


5 


X 22 = 0, det(x Q/3 ) 7^ 0, rc(a:) = k (j8 + a; 2 )" 4/3 


Yi, 3/?y 2 + Y4 


G 


X 22 7^ 0, det(x a/3 ) = 0, 

k(x)= K0 ((3 + xi-( X 12 /x 22 )x 2 r 4 


+ y 3 , 
(x 12 /x 22 )n + r 2 



7 


X 22 = 0, det(x Q/3 ) 7^ 0, 

k(x) = k (f3 + 2x 1 -( X 11 /x 12 )x 2 y 4 


m + 2Y 4 , 

(X n /X 12 )F1 + 2Y 2 


8 


X 22 det(x a/3 ) 7^ 0, k(x) = const 


Y u Y 2 


9 


X 22 det(x Q/3 ) = 0, «(ar) = const ^ 


Yu Y 2 


10 


X 22 7^ 0, det(x Q/3 ) = 0, k(x) = 


Yu Y 2 , Y 3 


11 


X 22 = 0, det(x a/3 ) 7^ 0, k(x) = 


Yu Y 2 , Y 4 



Having completely solved the group-classification problem, our next step is to identify 
the variational symmetries of those equations of form (|58l ) which have been found to admit 
larger symmetry groups. For this purpose, we are to apply the condition (|2^) to the linear 
combinations of (|14| ) and the vector fields presented in Table 1, the respective functions A 
being given by (|65[). Thus, we found that all vector fields quoted under numbers 1, 3, 5, 7, 
8, 9 and 11 generate variational symmetries as well. In case # 2 the variational symmetries 
are associated with + 2{ X 12 / X 22 )P 2 ] *i + 2/3 2 Y 2 + Y 3 + (1/2)X . Similarly, in case # 4 
the variational symmetries are generated by Y\ and 2f3Y 2 + Y 3 + (l/2)Xo, in case # 6 - by 
(x 12 /x 22 )Fi + Y 2 and 0Yi +Y 3 + (1/2)X , and in case # 10 - by Y x , Y 2 and Y 3 + {1/2)X . 

Once the variational symmetries of the equations (58) are identified, we can derive the 
corresponding conservation laws. The conserved currents of the conservation laws for the 
equations given in Table 1 are computed using (pfl) , the above notes concerning the corre- 
sponding variational symmetries being taken into account. The obtained conservation laws 
are listed in Table 2 (in the same order as in Table 1) in terms of the differential functions: 

B(i) = -\[l{^w 1 w 111 ~w 2 ll )+x 11 w 2 1 ~x 22 wl + nw 2 ]-\{x 2 ^ww ll ), 2 , 



(l) — _ 2 L H 1 111 — 11^ ~ A l <v 2 -rrvw\ — -\ 

B 2 1} = -X 2M WiW M + ^(x 2m ww m ),i 



B\ 2) = ~x 1 '"w 2 w t _ l . + j(wuWi2 ~ w 2 wui) - ^(ywiwu - x 1 * 1 ^), 2 , 
B 2 {2) = -^{Wn+X^wl-x^wl + nw 2 ) +^(jw 1 w 11 ~x ltl ww fl ), 1 , 
B^ 3) = (x 1 + ^x 2 ^ Bfc + 2x 2 B? 2) + x^ww^ + ]p5 la {ww xll - w lWll ) 

V 11 \ 



Bfa = [^ + ^)B^+3x^ 2) 



[x afi w Wf , + S la ( x 11 ww 1 + 2 x 12 ww 2 - jwtWn)]. 



Table 2. Conservation laws for equations of form rt5 



M. 

W 




Conservation laws 






1 


V 1 


)+ /? 2 i%]=0 






2 


DJiB 1 + 


2(X 12 /X 22 )/? 2 ) Bfa + 2(3 2 Bf 2 


+ 


B r3)]=o 


n 

6 


D a [(P + 


3(X U /X 12 )/? 2 ) Bf x) + 60 2 Bf 2 


+ 


2B ( * 4) ] = 


4 


D a B? x \ = 


0, D a [2pBfc + Bfa] =0 






5 


D a B"^ = 


0, D a [3(3B^ 2) +Bfo] = 






6 






+ Sf 2) ]=0 


7 


D a [8B^ 


+ 25« 4) ]=0, D a [( x u /X 12 ) 


Sf 1) +2Bf 2) ]=0 


8 




0, D a Bf 2) = 






9 


DaB"^ — 


0, D a Bf 2) = 






10 


D a B^ — 


0, D a B" 2 j — 0, D a B" 3 j = 







11 


D a B%s — 


0, D a B" 2 ^ — 0, D a B^ = 








According to the general results of Section 3, each equation (|5§|) admits conservation laws 
with characteristics Q = u(x), where u(x) is any smooth solution of the equation considered. 
These conservation laws are of the form d3fl) , that is 



D«Pfa = 0, 

the corresponding conserved currents P,", being given by the expression (^5|). Here, on 
account of (p9h, (B5h simplifies and reads 



P (u) = X"* 1 !^ - u ni w ) + ^ lQ 7( uu; iii + ".11 wi - u^u io - u,i ton). (80) 

Let us now specialize to the differential equation 

EJwuix + mu>22 = 0, (81) 

governing the dynamics of a classic homogeneous Bernoulli-Euler beam. Here EJ is the 
bending rigidity of the beam and m is the mass of the beam per unit length. According to 
the above analysis, (|8l| ) admits the following six linearly independent infinitesimal variational 
symmetries: 

Y u Y 2 , Y 3 + l-Xo, Y B = ^-, Y 6 =x 1 -^-, Y 7 = x 2 ^-, (82) 
2 ow ow ow 

where Y5, Yq and Y7 are vector fields of the type X u = u(x)d/dw corresponding to the 
solutions u = 1, u = x 1 and u — x 2 of (|8l|), respectively. Here, the independent variables 
x 1 and x 2 are the spatial variable along the rod axis and the time, respectively, so that the 
conservation laws admitted by the smooth solutions of equation (81) may be written in the 
more familiar form 

d<P dP _ Q 
dx 2 dx 1 

where W and P denote the density and flux of the conservation law, respectively. The 
densities and fluxes of the conservation laws for (^l|) associated with the vector fields ( |S2"|) 
together with their physical interpretation are presented in Table 3. 



Tabic 3. Conservation laws for Bcrnoulli-Eulcr beams 



Generators 


Conservation laws 


space translations 
Yi 


wave momentum 

$m = mwiW2 

P(i) = U/ ^)[EJ \£w\W\\\ — wfi) — rmofj 


time translations 
Y 2 


energy 

*/ 2 i = (1/2) {EJw 2 , + mw 2 ) 


P(2) = EJ{wuW\2 - W 2 Wlll) 


scaling 


*( 3 ) = a; 1 *^) + 2x 2 * (2 ) + TYIWW2 


Y 3 + (l/2)X 


P (3) = z 1 ^!) + 2x 2 P (2) + (1/2)^7(^1X1 - tuiiwn) 


Y 5 


linear momentum 


* (5 ) = mw 2 , P( 5 ) = EJwxu 


Y 6 


similar to Eshelby energy-momentum tensor 


W( 6 ) = X 1 miV2 1 P( 6 ) = PJ(x 1 wm - Wu) 


Galilean boost 


center-of-mass theorem 


Y 7 


^(7) = m{x 2 W2 — w), P(7) = EJx 2 W\\\ 



Conservation laws in the dynamics of rods are considered in many papers (see e.g. 
Antman, 1984; Kienzler, 1986; Chien et al., 1993; Maddocks and Dichmann, 1994; Tabarrok 
et al., 1994; Djondjorov, 1995). However, the particular form of the differential equations 
examined in the present study allows comparison with the results reported by Chien et al. 
(1993) and by Maddocks and Dichmann (1994) only. 

Chien et al. (1993) derive conservation laws for the statics and dynamics of rods em- 
ploying a technique called by the authors Neutral Action (NA) method. The conservation 
laws for rod equations established in this Section could be compared to their ones only for 
the differential equation ( |8l| ) which coincides with the equation 

Pwim + 2P iWin + P iiWn + Hw 2 2 = 0, 

considered in Chien et al. (1993) when B = EJ and H = m. The comparison shows 
that the conserved currents of the conservation laws for (|si|) with characteristics other than 
Q = u(x) obtained by Chien et al. (1993) coincide with ours presented in Table 3. As for 
the conserved currents of the conservation laws fo r (|8l| ) with characteristics Q — u(x), where 
u(x) is any solution of (|si|), our general formula (|8C|) implies 

P>"j = 8 2a m{uw2 — u,2 w) + 5 la E J {uwni + u, n wi — it,m w — u,i Wn). 

Only a part of these conservation laws are identified and presented in (Chien et al., 1993), 
namely those associated with the solutions to (|8l]) of the form 

u(x) = dix 1 ) 3 + C 2 {x 1 ) 2 + C3X 1 + C 4 + C 5 x 2 , Ci= const (i = 1, . . . , 5) , 

while, in fact, equation ( |8l| ) has an infinite-dimensional space of solutions. 

Five conservation laws in the dynamics of rods are reported in (Maddocks and Dich- 
mann, 1994) within a general nonlinear direct theory. The restricted version of this theory 
describing small planar bending of an uniform inextcnsiblc unshcarable isotropic elastic rod 
with a linear constitutive law, the rotatory inertia of the rod cross section being neglected, 
is exactly the classic Bernoulli-Euler theory for homogeneous beams whose governing equa- 
tion is (|Sl|). Rewriting the conservation laws in (Maddocks and Dichmann, 1994) taking 
into account the aforementioned restrictions we observe that: (i) the conservation law for 
the total angular momentum (formula 2.14 in Maddocks and Dichmann, 1994) degenerates 
to the well known basic relation of Bernoulli-Euler theory Q — dM/dx 1 (here Q and M 
denote shear force and bending moment, respectively, see Washizu, 1982); (2) the density 



and flux of the conservation law associated with material isotropy (formula 4.5 in Maddocks 
and Dichmann, 1994) vanish identically; (3) the conservation law corresponding to material 
homogeneity (formula 3.2 in Maddocks and Dichmann, 1994) reduces to conservation of the 
wave momentum (see Table 3); (4) the expressions for the densities and fluxes of energy 
(formula 2.19 in Maddocks and Dichmann, 1994) and linear momentum (formula 2.12 in 
Maddocks and Dichmann, 1994) conservation laws coincide with the respective ones pre- 
sented in Table 3. The set of conservation laws with characteristics Q = u(x), where u(x) is 
any solution of (|si|), as well as the conservation law associated with the variational scaling 
symmetry Y3 + (l/2)Xo (see Table 3) have no analogues in (Maddocks and Dichmann, 1994). 

Three interesting kinds of group-invariant solutions to certain equations of the class ( p8[ ) 
arc identified below. First of them corresponds to vector fields cYi =F Y2, where c = const. 
These group-invariant solutions are travelling waves 

w = U(s), s = x ± ex 2 , 

admissible only for equations ( ]58| ) with k^x 1 , x 2 ) = /(s). The reduced equations determining 
such group-invariant solutions are 

7^ + (X 11 ± 2x 12 c + X 22 c 2 )^ + f(s)U = 0. 
The second one corresponds to the vector field Y3 and is of the form 

w = U(s), s = x 1 {x 2 )- 1 ' 2 ~^ I {x 2 ) 1 / 2 . 

X 

The vector field F3 is admitted only if ^(x 1 ,^ 2 ) = (a; 2 ) _2 /(s) (see cases # 2, 4, 6 and 10 in 
Table 1). The reduced equations for these invariant solutions are 

d 4 U 99 9 d 2 U 99 dU 

4 ^ + * s V +3 ^ +4/ ^ = - 

The third kind of group-invariant solutions corresponds to the vector field Y^: 

w = u(s), s = 2xVr 1/3 -C(* 2 ) 2/3 . 

X 12 

The vector field Y4 is admitted only if k(x x ,x 2 ) = (a; 2 )~ 4 / 3 /(s) (see cases # 3, 5, 7 and 11 
in Table 1). The reduced equations for the invariant solution under consideration are 

Obviously, the latter two kinds of group-invariant solutions could be reduced to self-similar 
solutions if x 12 = or x 11 = 0, respectively. 



6. CONCLUDING REMARKS 

In this paper, Lie transformation group methods have been applied to the class of partial 
differential equations ([]]) . This class is of interest for structural mechanics since the governing 
equations of various classical plate and rod theories belong to it; the examples given in 
Section 2 illustrate this fact. In the context of structural mechanics, the results of the group 
analysis of equations ([!]) give a number of attractive possibilities. Here, the established 
point Lie symmetries of (Q) are used to construct group-invariant solution to the governing 
equations of several plate and rod models, to derive conservation laws revealing important 
features of such models and to find transformations simplifying the differential structure of 
equations associated with particular plate problems. 



First of all, the well known computational procedure for finding the most general point 
Lie symmetry group has been applied to the foregoing class of equations. As a result, the 



system of equations (16) - (^0|) is derived determining the equations of the type considered 
that admit a larger group together with the generators of this group; naturally, all equations 
of this class being linear and homogeneous admit the point Lie groups generated by the vector 
fields (|l4|) . The system ( |l6| ) - ( ]20j ) allows the associated group-classification problem to be 
stated and examined. 

In Section 4, this problem is solved for the plate equations ( |38| ) in terms of their invariants 



sm, S( 2 ) and srg\ defined by (47). The equations of form (38) with sm = sw\ = S(3) = are 
found to admit the largest symmetry groups. It is noteworthy that each such equation with 
variable coefficients can be transformed, using a suitable change of variables, to an equation 
with constant coefficients belonging to the same class. An example of such a transformation 
is given at the end of Section 4 where, in addition, a class of group-invariant solutions to 
the equation considered is presented. 

The group-classification problem for the rod equations ( [58|) is entirely solved in Section 5. 
All equations of that kind admitting point Lie symmetry groups, in addition to the ones 
generated by |l4|), are determined and presented in Table 1 together with the generators 
of the respective groups. The largest symmetry groups are admitted by the equations of 
form ( p8| ) whose coefficients are such that x 22 det(x ct ' 3 ) = 0, k(x) = 0. The most interesting 
group-invariant solutions for equations ( |58| ) are identified and the corresponding reduced 
equations are presented at the end of Section 5. 

Once the "ordinary" point Lie symmetries of an equation of form (|l|) are determined, one 
can easily find, using the general criterion ( p^ ) , which of them are variational symmetries of 
this equation. Then, ([35|) and (|37]) provide explicit expressions for the conserved currents of 
the conservation laws associated through Noether's theorem with the established variational 
symmetries. These expressions will involve derivatives of the dependent variable of lowest 
possible order, which is important in view of their application in structural mechanics. The 



is given explicitly by formula (p3|). 
27|) of the general criterion (|26|), that 



reciprocity relation valid for each equation of form 

In Section 4, it is shown, using the consequence 
each point Lie symmetry of a plate equation of form ( 38ft generated by a vector field of form 
( [46| ) is variational symmetry of this equation. Therefor, each such symmetry gives rise to a 
conservation law with characteristic Q = (l/2)^w — and conserved current given by 

( j37l ) admitted by the smooth solutions of the respective equation. 

The conservation laws for the rod equations listed in Table 1 are given in Table 2. 
Inspecting these results one can see that the equations for unsupported rods and rods on 
Winkler foundations admit two independent conservation laws associated with the wave 
momentum (D a B^ = 0) and energy (D a B^ = 0). Equations @ and ([!(]) governing 
the stability of unsupported axially compressed beams and fluid conveying pipes belong 
to this class. Rod equations with k(x) = and det(x Q/3 ) = admit a supplementary 
conservation law D a B^ = associated with the infinitesimal scaling symmetry Y3. Such 
an equation is ( |8l| ) governing the vibration of the classic Bernoulli-Euler beam. Table 3 
contains several physically important conservation laws for this equation. A comparison 
between the conservation laws derived here for equation (^l|) and the relevant results in 
(Chien et al., 1993) and (Maddocks and Dichmann, 1994) is presented in Section 5. 

Acknowledgements — This research was supported by Contract MM 517/1995 with the NSF, 
Bulgaria. 



REFERENCES 

Antman, S.S., 1984. The theory of rods. In: Mechanics of solids, Vol. II. Springer- Verlag, 
Berlin, pp. 641-703. 

Bluman, G.W., Kumei, S., 1989. Symmetries and Differential Equations. Springer- Verlag, 
New York. 



Chicn, N., Honcin, T., Herrmann, G., 1993. Conservation laws for nonhomogcneous 
Bcrnoulli-Euler beams. International Journal of Solids and Structures 30 (23), 3321- 
3335. 

Djondjorov, P., 1995. Invariant properties of Timoshenko beams equations. International 
Journal of Engineering Science 33 (14), 2103-2114. 

Gregory, R.W., Paidoussis, M.P., 1966. Unstable oscillation of tubular cantilevers convey- 
ing fluid. I: Theory. Proceedings of Royal Society London A-293, 512-527. 

Ibragimov, N.H., 1985. Transformation Groups Applied to Mathematical Physics. Reidel, 
Boston. 

Kienzler, R., 1986. On existence and completeness of conservation laws associated with 
elementary beam theory. International Journal of Solids and Structures 22 (7), 789- 
796. 

Maddocks, J., Dichmann, D., 1994. Conservation laws in the dynamics of rods. Journal of 
Elasticity 34, 83-96. 

Olver, P.J., 1993. Applications of Lie Groups to Differential Equations, Second Edition, 
Graduate Texts in Mathematics, Vol. 107. Springer- Verlag, New York. 

Olver, P.J., 1995. Equivalence, Invariance and Symmetry. Cambridge University Press, 
Cambridge. 

Ovsiannikov, L.V., 1972. Group properties of the equations of mechanics. In: Mechanics of 
Continuous Media and Relevant Problems of Analysis. Nauka, Moscow, pp. 381-393. 

Ovsiannikov, L.V., 1982. Group Analysis of Differential Equations. Academic Press, New 
York. 

Smith, T.E., Herrmann, G., 1972. Stability of a beam on an elastic foundation subjected 
to follower forces. Journal of Sound and Vibration 39, 628-629. 

Tabarrok, B., Tezer, C, Stylianou, M., 1994. A note on conservation principles in classical 
mechanics. Acta Mechanica 107, 137-152. 

Vassilev, V., 1988. Group properties of a class of fourth-order partial differential equations. 

Annals of the University of Sofia 82 (Part II - Mechanics), 163-178. 
Vassilev, V., 1991. Group analysis of a class of equations of the plate and shell theory. Ph.D. 

thesis, Institute of Mechanics and Biomechanics, Bulgarian Academy of Sciences, Sofia. 
Vassilev, V., 1997. Application of Lie groups to the theory of shells and rods. Nonlinear 

Analysis 30 (8), 4839-4848. 
Washizu, K., 1982. Variational Methods in Elasticity and Plasticity. Pergamon Press, 

Oxford. 



